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Orthomorphic Projection of an Ellipsoid upon a Sphere, 

By Thomas Craig, 

Johns Hopkins University, and U. S. Coast and Geodetic Survey. 

After the labors of such men as Gauss, Lagrange, and Lambert, upon the 
general theory of orthomorphic projection, or projection by similarity of infini- 
tesimal areas, it does not seem as though much of value or interest could be 
obtained by any further study of the subject, and indeed there is nothing more 
to be said upon the general theory, but an abundance of opportunities still exists 
for applying the results obtained by these princes of the realm of mathematics 
to the solution of particular cases of the problem. Boole, in the supplementary 
volume to his Differential Equations (a posthumous work, edited by Tod- 
hunter), gives a most elegant investigation of orthomorphic projection upon a 
plane ; the formula? arrived at are applicable to any surface which it may be 
desired to project upon a plane. One application is made by Boole to the case 
of an oblate spheroid such as the earth. I have given in another place an appli- 
cation of his formulas to the projection of the general ellipsoid upon a plane, 
and the results obtained there led me to attempt the projection of the ellipsoid 
upon a sphere. Gauss gives, among the examples illustrative of his general 
theory, the projection of an ellipsoid of revolution upon a sphere, and considers 
the so-solved problem to be a valuable addition to geodesy. At the present day, 
when it is at least considered possible that the earth may be a general ellipsoid, 
it may be that the formulas necessary for its projection upon a sphere will not 
be devoid of interest. 

Denote by R the radius of the sphere, then its equation may be given as 

£2 + ^ + £2 =jR2 . (1) 

if U and V are the usual spherical co-ordinates, we have for g, rj, £, the values 

£ = R cos U sin V, 

v = R sin U sin V, (2) 

£ = R cos V. 
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Denote by a, h, c, the serai-axes of the ellipsoid, then this surface is given by 

o o o 

— * + y + -j — 1 ; (3) 



if X t and X 2 denote the variable parameters belonging to the two hyperboloids 
confocal to the given ellipsoid, we have for these surfaces the equations 



a 2 + \ ^ b 2 + A, ^ c 2 + \ ~ 1 > 

"T" 7,2 i \ ~T „2 , \ ij 



(4) 



(7) 



a 2 + A 2 ^ b 2 + A 2 ^ c 2 + A 2 

and then, as is well known, the co-ordinates x, y, z, are given by 

, _ a 2 (a 2 + A,) (a 2 + A,) 

* ( a 2 _ #>) (a 2 - c 2 ) 

2 _ 5 2 (5 2 + A,) (b 2 + A 2 ) 

# — (^ _ e ") (J2 _ a 2) i,Oj 

2 _ c 2 (c 2 + A.) (c 2 + A 2 ) 
~ (c 2 - a 2 ) (c 2 - 6 2 ) 

and the element of length on the ellipsoid by 

" 4 \(a 2 + \) (b 2 + A x ) (c 2 + A,) ^ (a 2 + A 2 ) (b 2 + A 2 ) (c 2 + A 2 ) J ' W 

Write for brevity 

L x = V(a' + A 1 )(ft , + A 1 )(o ,, -|-.A 1 ), 

Z 2 = V(« 2 + A,) (6 2 + A 2 ) (c 2 + A,) ; 
the Gaussian equation 

a = 

leads then to 

A A ' 

the differential equations of the problem. 

For the element of length on the sphere equations (2) give 

& = dS 2 = R 2 sin 2 VdU 2 + i^F 2 . (9) 

Q} — leads to the differential equations 

and from this follows by integration 

27" =b i log cot 1 V = const. (10) 

If we find one integral of (8), in the form 

P 4- iQ = const., 
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the problem will be solved by equating U to the real, and i log cot \ V to the 
imaginary part of 

f{P + iQY; 

f being an arbitrary functional symbol and P + iQ being the argument of the 
function, all operations having to be performed upon this complex-quantity as a 
whole and not upon its components. 

The first thing to be done is, of course, to so transform (8) that an integration 
may be performed. To this end, observe that the parameters X x and X 2 are 

limited by the relations 

- c 2 > X, > - b\ 

— 6 2 > X 2 > - a\ 

so that we can express Xj in terms of a new variable 6 as 

b 2 (a 2 - g 2 ) cos 2 + c 2 (a 2 - b 2 ) sin 2 d . 



Xi = 



or again, writing 



{a 2 — c 2 ) cos 2 6 + {a 2 — b 2 ) sin 2 



a 2 - & . iQ 2 

^F^7 2 tan e = °>> 



Xi = 



b 2 + G 2 



1 + c, 2 



Prom these we have the relations 



d 2 +X 1 



(a 2 - b 2 ) + (a 2 - c 2 ) 






1 + 

(6 2 - e 2 ) to 2 

(b 2 - c 2 ) 



I 



A, dX t 



K = - 



b 2 + c 2 

1 + w s 



2~ ? 



- 2 - C 2 ) 2 oxfa 
"*1 — (1 + w 2y2 



These substitutions made in P give after simple reductions 

clta 



p 



f: 



Vb 2 + cV 



V(a 2 - b 2 ) + (a 2 - c 2 ) <o 2 1 + "> 2 

For the transformation back to the variable 6 'observe that we have 

b 2 (a 2 - e 2 ) - a 2 (b 2 - c 2 ) sin 2 6 1 



J2 + ^2 = 



(a. 2 — e 2 ) cos 2 



do = \/ a \~ h \ sec 2 d0, 

a 2 r- c 2 

(* 2 - o 2 ) cos 2 6 



1. + cu 2 = 



(a? - g 2 ) - (b 2 - c 2 ) sin 2 61 ' 



V(a» -**) + («* - c 2 ) co 2 = ^1_J1 ; 
v cos 6" 



(11) 

(12) 
(13) 



(14) 



(15) 



(16) 



Craig : Orthomorphic Projection of an Ellipsoid upon a Sphere. 117 

multiplication of the last three of these gives 

2 V« 2 — e' 2 cos Ode 



(V - c 2 ) - (6 2 - c 2 ) sin 2 6 
The result of substituting in P these values is 



/ 


L 1 6 2 (« 2 - C 2 ) Sm ^J 




6 2 — e 2 
1 2 2 sin 2 # 



p = /_jl ? / '^ :;yjo"° -i J * d7) 

The quantity 



a 2_ c a< 1 

occurs in both numerator and denominator of the differential expression ; in the 

a 2 
numerator occurs also the factor — ; for a 2 = b 2 , 

ft 2 b 2 -C 2 _ 

b 2 ' a 1 - c 2 ~~ 1 > 

for <? = £ 2 , this is = ; so that the- value of 



V 1 a 2 -c 2 



lies always between the limits and 1. 
Write then 



52 a 2 yi _ 2 yi _ c 2 ■ 

5 = sin 2 a, ^ • ^372 = A 2 , ^— 2 = P sin 2 a , (18) 



a 
and P takes the form 

\em 



_ Bin a Aq /• Afl dg 

cos a J 1 — l£ l sin 2 a sia 2 # ^> ' 



since 



Aa = l/*3 and cos a = V^-E* 






If we call e the eccentricity of the section of the ellipsoid made by the plane 
xy, it is clear that 

a = cos- 1 c. (20) 

The above value for P can also be written in the form 

P = tan a Aa J ^ - F sin a cos a Aa J (1 _ ^ sin 2 g sin 2 e) A0 ' (21) 
Introducing elliptic functions by means of the equations 

C d& fd* 

t= J 55' « = J Aa 
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r k 2 sn a en a dn a sn 2 t . dt 

(21) becomes P - tn a dn a .t - J ± - k * &n * a sn * t ' (22) 

The quantity under the integral sign is (Cayley's Elliptic Functions, page 15) 
the elliptic integral of the third kind, or n (t , a) ; but introducing Jacobi's 
functions Z and © by the relation 

n(*,a) = *Z» + *log§g=^, (23) 

we are enabled to write immediately 

P = [tn a dn a - Za\t-\ log ||~j , (24) 

or, writing log e 2 [tn ° dn a ~ m l = log e 2 *", 

'-*l»gj^-'-- (25) 

For the value of w observe (Cayley's Elliptic Functions, page 157) that 

H (« + K) = © («) |/| en «, 



and further that 



<Z snw dn m . 

^ log en u = cnu =tnudnu 






"We have then 



^ log H (a + -ff') = tn « dn m — Z (w) 



The functions © and H expressed in terms of the ^-functions are 

(~V*) = 1 — 2 ^ cos 2 r + 2 ? 4 cos 4 r — 2 / cos 6 t + . . . 

H ( — ^ J = 2 f/# { sin r — # 2 sin 3 r + # 6 sin 5 7 — j 12 sin 7 r + . . . 

or, writing — f = t, these may be expressed as 

e (0 = e (^T *0 = 1 + 2 S (- iy ^ cos 2/^ 

H 00 = H (v f ) = 2 ^ 2 (- I)''" 1 qM-v sin (2/- 1) f, 

and consequently 

© ft + a) _ 1 + 2 S (- I/ y* cos 2 j (g + ftp 
®(t-a) ~ 1 + 2 2 (- iy y* cos 2 j (*' - a') 
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and H (a + K) = 2 .fq t g*U-i) cos (2/ - 1) a', 

where a — 2k • 

Since d _ ■* d 

da 2 K da 
we have for w the value 

w = ~ da lo S H (« + K ) = ~ YK da lo S H (a + K) 

__ _ _tt_ 2 (2j - 1) g"-'- 1 ' sin (2j - 1) a' 
~~2if 2 ^'< J - 1) cos (2j -1) a' 

The complete value of P is now 

p-il C ® g - a ) ** 2 (2j - 1) gJQ-i) sin (2j- 1) a / ) 

■^ — 2 l0 S |© (« + a) ' ex P- if 5 j/iy-i) C os (2,/ - 1) a' ) ' 

exp. « standing for e w . 

If, instead of the quantities t and a, their complements are introduced as 

K — t = 4 , j5T — « = «! 
and also 

am (ifj) = #j , am (a^ = a x , 

then the integral expression for P becomes 

_ sin aAa f AOdO 

cos a J 1 — k 2 



_ cos 
sin 

Now 



in 04 J (1 - /c 2 sin' 2 a x sin 2 0,) A(9 X ' 

«^ _ ^1 pi + & sin 2 a, sin 2 0,-| , 

A6>! (1 - A 2 sin 2 «j sin 2 (9J ~ A^ |_1 - W 2 sin 2 a, sin 2 (9, J ' 

and therefore, by changing the order of integration and dropping a useless 
constant, 

__ cos q.Aa, Cdd x r dd x 

P — g ; n ~ — / -TfiT + " cos a x sin aiAa, I -rzTTi 72 • 2 --"O^ 

sin a x J A(/j l 1 1 J A0J (1 — & sin 2 a : snr #J 

_ dlogH(a l ) _ (R) ft, + a t ) 

rfo, ? i 2 l0 g © ( tj _ a ,) ' 

If we introduce an imaginary argument and the complementary modulus, this 
expression takes another form rather more convenient for computation. We 
have (Cayley's Elliptic Functions, page 151) 



%{u) = \/ ] ^ **!£ A- e (*«,#), 
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or again © (u) = j/|i e~*KK' H (K' ± iu , //) 

and H («) =—y±e*x*U (iu , Jc') ; 

now, denning w l as 

it is obvious that 

w = ~ l 2KK> + w o= ^ % Hfo)- 
Now, introducing the ^-functions, 

H (ifflj. , k') = 2 ^' {sin w^ — j? sin 3 «'«{ + yj sin 5 ia[ . . .} 
and therefore, introducing exponentials, 

_ j^ % (_ iy-i (2j - 1) gl "^-" [e^- 1 '".' ■+. g-w-W] 

where «{ = ^^ . The function P has already been. determined as 



2 if 



P - < - 1 log 0(<i _ Oi) > 
and can now be written as 

■^ W/l 2 io g H (J? - i (i L - a,) , k') 
The logarithm of this ratio of conjugate H-functions is given by 

where t[ = ^-gj • The final form of P is then 

_ jr^_ 2 (- l) J - 1 (2j - 1) y^'" 1 ' [ e «-iw + e -<« -««i'] 
~ 2f 2 (- l) J '- 1 (2y - 1) y''^- 1 ' [ e w -•)«>' _. e -(2j -««»'] 



l 



log; 



Jj/U-1)MJ-1|(1i'+«i']^ g-VJ -1)W +ai')l 



2 lv, o jA/U-lir^-ilW-m') e -(2j'-D((i' — 0,'n 

By interchanging X x into A 2 the integral expression for P becomes 

The quantities It , a do not depend on either \ x or X 2 , and in consequence are 
unaltered by this change; the same is, of course, true of the constants a,d,w, 
which are functions of Tc , a and other constants. The only quantity which can 
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vary is t, and this, on account of the prescribed limits of X 2 , will become a pure 

imaginary, say 

t = i(E' + r); 

then = am i (JT + t) ; 

now stfi{K' + T)= w ^r Tr 



and 
but 



i sn (t , 7cQ 
sn(fT,*)- cn(TjA0 itn(r,*)j 

1 _ (a 2 + A,) Z, 2 _ 1 . 

k l sn' 2 « ~ a 3 ' (b 1 + A 2 ) ~ ¥ sn 2 i (JT + t) ' 

therefore tn 2 (t , //) = — sn 2 («' t) = ft2 /^ _ A • 

Writing i/> = am (t , //) , 

-,. . . ft 2 (S 2 + A,) . (ft 2 + A s ) &» 

we derive cos 2 ^ = (a ._p )A , . «m 2 ^ = (o *_j. )Af > 

and,. for the complementary modulus &', 

_ c«( a '-y) 

A — 6' 2 (a 2 -c 2 )' 
from which follow #' 2 sin 2 ii = 7-5 |rr- 

and A (i/» , £') = /l - 7/ 2 sin 2 xjj = |/^+A> . 

(ft — ) A 2 

Now, since P = «rf + log j/® (*+_") 

and Q = i P> 

it follows at once that 

£ - 2* (K + t) - g log @( i K , + iT _ a) > 

but © (« + »£") = * e 4 ^ ( *' _ 2 iU) H («) ; 

therefore = w (JT + t) - | log |^j ^ , 

or, dropping the constant iv K ', 

— _f -1 2(-l)'-V u -"cos(2j-l) a' [e' 2 '- 1 ^ -e-o-D*] 
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where a = g-g., t = ^ • 

On transforming Q to the complementary argument a t it becomes 

£ - »r - 2- log H(j ^_ ai _ iT) 

1 , © (t + m, , &') 
= w lT - 2". log 0(T _ t - Oij ^ > 

or finally 

Q — w x r — tan f_ 2 (- Vf <f cos 2^ T ' [eW - e- 2 -""'] * 

We have said that it is only necessary, in order to solve the proposed prob- 
lem in the most general manner, to place U equal to the real and i log cot ^ V to 
the imaginary part of 

/(P + iG), 

where / is an arbitrary functional symbol ; this is the same thing as writing 

U + i log cot i V=f(P + iQ), 

U— » log cot | V=f{P — iQ). 
Suppose that we take 

f{v) = v, 
then follow 

U=P, 

F^COt- 1 ^; 

by the first of these relations all curves which depend only upon t are projected 
into curves depending only upon U, the longitude ; that is, all the curves de- 
pending upon t are projected into meridians ; but 



t 



~J AO' 



and 6 is a function of X t only, so that t is a function of \ 1} and the curves which 
depend upon t are the lines of curvature cut out of the ellipsoid by the hyper- 
boloid of two nappes; these lines of curvature are then projected upon the 
sphere in the meridians, and in like manner the equation 

V= cot" 1 e Q 

shows that the second system of lines of curvature is projected in the parallels 
of latitude. 

If the sphere be projected into the plane £r? by the relations 

(U+i log COt JF) = log (f + i-q) 

(U— i log cot |- V) = log (f — iy), 
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all of the parallels will be projected in concentric circles, and all of the meridians 
in straight lines passing through the centre of the circles (the stereographic pro- 
jection). Now, if we introduce polar co-ordinates {p , </>) by the formulas 



we have 

but 

and 

therefore 
but 



consequently 



$ = p cos (f> , 7) = p sin (j>, 

£2 _| 2 — „2 . 

£ -1-17 — p ; 

log (f + itj) + log (f — iyj) = log (f 2 + ?f) = log p 2 

i / ,*- i • \ i if • \ i cos <£ + i sin <jf> 

log (f + tr,) - log (f - .,) = log coa<ft _ iain<ft = 2 ^<£; 

F= cot -1 e* ; 
F=cor 1 ^, 



P = w t J ® (t + a) 
H V © (* - a) ' 

• _ /i _ ^ _i_ 2 i H (« + * r ) 

<k-Q-™t+ Ti bg R(a _; r) • 

We have in this case the lines of curvature on the ellipsoid arising from its inter- 
sections with the hyperboloid of two nappes projected on the plane in straight 
lines passing through the origin of co-ordinates ; the remaining system of lines 
of curvature being projected in concentric circles whose common centre is at the 
origin of co-ordinates. If the ellipsoid be one of revolution around the axis of x, 
then result 

b = c , k = o ,k' = 1 ,q = o , K= k , K' = o , 

8 = t — t' , a = a = a , w = tan a , 



■ / = T = /lc^ = 1 °g t ^H9o o + ^) 



and 



e T 



TT^ 7 ' = sm ^ 
For the position of the point x, y, z-, when projected on the sphere, we have 

U= tan a . 0, 

V= tan a . log tan J (90° + \f>) + tan -1 [cot a sin <£], 
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and for the corresponding co-ordinates on the stereographic projection 

p = e tana - e 
<f>=V. 

For the oblate ellipsoid we must use the second forms that we obtained for P 
and Q, and introduce the conditions 

h' = o , q' = o , K' = g > K— o,t = r = xfj , 

a[ = a t = log tan \ (90° + a*) ,f 1 = t i = log tan £ (90° + 6 l ) i 

^ tan J (90° + a x ) - cot A (90= + aj l0 ° tan 2 l yU ^ ^) 

1 tan & (90° + a,) tan | (90° + 6> 3 ) + cot j- (90° + a t ) cot | (90° + 6 X ) 

2 lo S cot A (90° + 0l ) tan £ (90° + 0J + tan £ (90° + oj cot £ (90° + ^) 

1 i x i /nno i a % i i sin ' 2 i Ci + a i) + cos2 i ( 6 i — tt i) 

- log tan I (90° + 6 t ) — I log — „ . ;' , , , . 5 . ; ' -{ 



§ 



sin 

-T— log tan i (90° + ft) — i lo<? T" ,. 

sm«i ° 2 v. x > 2 o 1 — sin a x sm X 

■A 



1 + sin a, sin ^j 



sin a, 



The projection on the sphere is determined by 



1 



TT , \ [sin A (90° + (9,)l sln "i Vl — sin a, sin 6, ) 
U=l0g i ±= i-j- \ \ , 

( [cos J- (90° + 5 1 )]Bina 1 |/l + gin ttj sin X ) 



F= cot -1 e sin «j ; 
that on the plane by 



P = 



[sin J (90°+e i )yn H |/l 



sin a, sm t 



[cos J- (90° + (9 x )]sin «, |/l + sin a, sin Q x 



4> = cot * e si » «i . 

For the complete solution of the problem it now remains only to give the values 
of x, y, s, the co-ordinates of any point on the surface, in terms of the same 
variables that have been employed in finding P and Q. We had 



2 __ a? (a* + xj ( ffi 2 + -x,) 

X ' (a 2 - V) (a 2 - c 2 ) ' 

, 5 2 (& 2 + A,) (5* + A 2 ) 
* (6 2 - c 3 ) (6 2 - a 2 ) ' 

„2 - c ° (? + X.) (fi» + A,) , 

g — ( c * _ a 5) ( c 2 _ p) ' 
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. 2 , __ (a 2 — 6 2 ) sin 3 a cos 2 t/f 
+ A2 — f. — J 



now 


.. | x _«?-V 




a ' 1 — Z> ' 




7 9 , x (a 2 — 6 2 ) sin 2 a sin 2 # 
* + Aj — ^ 




. (a 2 — c 2 ) sin 2 a cos 2 6 
c 2 +.K= 2 > i 


where 


2?, = 1- 



c 2 + X 2 = 



(a 2 -c 2 )sin 2 aA 2 0,F) 



A 



Z> 2 = sin 2 a + cos 2 a sin 2 *|», 

or, introducing the notation of elliptic functions, 

6 = am t , xf/ = am (r , //), 

Z?i = 1 — sn 2 a sn 2 1, 

D 2 = 1 — sn 2 « sn 2 (r , #') , 
and 

„2 , x _ »'-y . , x _ (^-^)sn 2 (x ; ^) 

^X^^sn^sn 2 *, JHX3= (^i^M), 

c 2 + X 1 = ^- ) sn 2 acn 2 ^, c 2 + X 2 = ^"^ "^ *»' (t ' *> . 

Observing now the relations 

ffi 2 -6 2 _ Q 2 -& 2 )5 2 _ a 2 A 2 a (» 2 -e 2 )o 2 _ « 2 7c' 2 

a 2 _ c3 — A a , //2 _ fl2 — F , J2 _ c 2 — A 2 a _ £2 ? 

and substituting in these dn a for A a, we are enabled to write x, y, s, in the 

forms 

x = G t . dn 2 a sn (t , #'), 

y = G x . dn 2 a sn 2 a sn en (r , #') , 
s = G t . 7/ sn 2 a en dn (t , A') , 



where 



fc= " 



dn a ^/(i _ A i sn a sn ^ ( Sn 2 a + cn 2 a sn !) ( T ? ;</) 5 

the a in the numerator of course denotes the semi-axis major of the ellipsoid. 
The angle a has been defined by 

from this 

b 2 = <# sin 2 a , 
and we also obtain quite readily 

„ a 2 k' 2 sin 2 a 

C ~ T2 5 

A a 
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the equation of the ellipsoid can thus be put in the form 

x 2 y 2 s 2 A 2 a _ 

^~ 2 "■" a 1 sin 2 a + a 2 k'* sin 2 a = - 1 ' 

or, transforming to elliptic functions, 



.?/ 2 



2 



* "^~ sn 2 a + A' 2 sn 2 a — a 

For the case of b = c it is easy to see that 

# = # 2 . sn (t , //), 
^ = (7 2 . sn 2 « sn ^ en (r , #') , 

= # 2 . sn 2 a en en (r , //), 
where 

G 2 = 



Vsn 2 (t , &') + sn 2 a en 2 (t , /c') ' 

For the case of the oblate ellipsoid it is necessary first to transform to the com- 
plementary arguments. This is done by means of the relations 



a 2 — c 2 



sin 8 aj = - , cos 2 a 1 = ^r, A 2 a t = ^ , 

sm ?i A^-c 2 )' 008 *i - A, (c 2 - 6 2 ) ' A *i ~ A, ( C 2 - a 2 ) ' 

The general ellipsoid is given in terms of the new elliptic function a x by the 
equation 

P dn 2 a s? 

x ~ en 2 a t ^ en 2 a 2 — a > 

and the co-ordinates jc, #, 2, of a point on its surface by 

x = G . dn 2 % dn ^ sn (r , k') , 
y = G . en 2 «i en t x en (t , #') , 
z = G . en 2 «! dn 2 a sn t x dn (t , £') , 



where 



G — 



a 



dna i V | (1 - A 2 sn 2 a, sn 2 f,) (en 2 a x + k rl sn 2 a, sn 2 (t , A') \ ' 

The oblate ellipsoid has b = a, and the equation of the surface becomes 

^ 2 + ?/ 2 H 2 — = a 2 , 

J en 2 «.j ' 

the co-ordinates x, y, z, being given by 

x = G 3 . en t x sn (t , Jc) , 
y=G s . en ^ en (r , £') , 
z = G 3 . en 2 «! sn 2 ^ , 
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where „ * t 

Vl — sn 2 o.j sn 2 t x 

The angle t/» is here the longitude and 6 1 the excentric anomaly of the meridians. 
The rectangular co-ordinates g, rj, £, can of course be given in terms of (t, a, r) 
or (t x , «i , r) by means of the relations 

f = R cos U sin V, 
7) = R sin U sin F, 
£ •= i? cos F; 
in terms of P and Q these become 

. R cos P 



7J — 



Ve 2<3 + 1 
PsinP 

R<? Q 



